Temperature and magnetic field dependent longitudinal spin relaxation in 
nitrogen-vacancy ensembles in diamond 
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We present an experimental study of the longitudinal electron-spin relaxation time (Ti) of nega- 
tively charged nitrogen-vacancy (NV) ensembles in diamond. Ti was studied as a function of tem- 
perature from 5 to 475 K and magnetic field from to 630 G for several samples with various NV 
and nitrogen concentrations. Our studies reveal three processes responsible for Ti relaxation. Above 
room temperature, a two-phonon Raman process dominates, and below, we observe an Orbach-type 
process with an activation energy, 73(4) meV, which closely matches the local vibrational modes of 
the NV center. At yet lower temperatures, sample dependent cross relaxation processes dominate, 
resulting in temperature independent values of Ti, from ms to minutes. The value of Ti in this limit 
depends sensitively on magnetic field and can be tuned by more than an order of magnitude. 

PACS numbers: 76.60.Es, 61.72.jn, 81.05.Uw 
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Negatively charged nitrogen- vacancy (NV) color cen- 
ters in diamond have been a focus of recent research 
due to their promise as fluorescent markers for biologi- 
cal systems [1 jSL qubits that can be optically initialized 
and readout J 1| , and sensors of magnetic and 
electric [lH fields. The NV center is uniquely suited for 
these applications due to its atom-scale spatial resolu- 
tion and exceptional optical and spin properties over a 
wide range of operating temperatures. 

The spin phase coherence time (T2) is a critical figure 
of merit for these emerging quantum-based applications. 
For example, for ensemble magnetometry, the sensitiv- 
ity scales as {NT2)~^^^ , where N is the number of spins. 



Therefore, optimization of ensemble sensors |ld. Il2|-|l4| 
calls for high-density samples with long phase coherence 
times. In quantum computing, T2 constrains the min- 
imum gate operation time and limits t he p erformance 
of quantum error correction protocols 15, IT6| . Typi- 
cally T2 is limited by magnetic impurities in the local 
environment which serve as sources of decoherence [l7| . 
However, it is possible to reduce the effect of the spin 
bath using dynamic- decoupling techniques [l3-[2lj. The 
phase coherence time is then limited by stochastic pro- 
cesses, such as phonon interactions, which cause irre- 
versible changes in axial spin projection. Such processes 
are known as longitudinal relaxation, with a character- 
istic timescale, Ti. 

The temperature dependence of Ti has been experi- 
mentally studied previously at high magnetic fields us- 
ing electron paramagnetic resonance j22l424| . In Ref. 
|22| . Ti was measured in the range 100-500 K, and 
the results were described by a model including Orbach 
and Raman phonon processes. Other results |23| over 
a somewhat lower temperature range could not be ac- 
curately described by this model. In Ref. Ti was 
measured at T « 2 K for two different samples, and the 



results differed by more than an order of magnitude, 
casting further doubt that Ti can be explained purely 
by phonon processes. 

In this Letter, we report an investigation of longitu- 
dinal spin relaxation Ti of NV ensembles as a function 
of sample impurity content, temperature, and magnetic 
field. We identify the principal relaxation mechanisms 
as interactions with local and lattice phonons, as well 
as cross relaxation with nearby spin impurities. We de- 
velop a universal model for Ti relaxation of NV centers 
which agrees well with experiments on samples across 
a wide range of impurity concentrations and tempera- 
tures. 



TABLE I. Sample characteristics (see Ref. |12 | for details). 
The diamonds were irradiated with 3 MeV electrons and an- 
nealed at 1050 °C for two hours. [NV] is the concentration 
of negatively charged centers and its estimate is accurate 
to within a factor of two. [N] is the concentration of para- 
magnetic substitutional nitrogen. HPHT - high-pressure 
high-temperature. CVD - chemical vapor deposition. 



Sample 


Synthesis 


Radiation dose 


[N] 


[NV] 


No. 




(cm-2) 


(ppm) 


(ppm) 


S2 


HPHT 


1 X lO'"* 


40-60 


16 


88 


HPHT 


4 X 10" 


10-30 


0.3 


S3 


CVD 


4 X 10" 


0.1 -4 


0.01 



The diamond samples used in this work are listed in 
Table HI These samples represent a wide range of nitro- 
gen (N) and NV concentrations. Measurements were 
performed using a confocal-microscopy apparatus de- 
scribed previously (lol . [2^ . An aspheric lens with focal 
length f = 6 mm and numerical aperture NA = 0.55 
was used to focus 532 nm laser light on a diamond sam- 
ple mounted in a flow-through liquid- helium cryostat. 
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FIG. 1. (a) Level structure of the NV center showing spin- 
triplet ground and excited '^E states, as well as the sin- 
glet system involved in the optical spin-polarization mech- 
anism, (b) Optical and microwave pulse sequence used for 
7i measurements, (c) Example of Ti measurement. The 
difference in level of fluorescence for two measurements with 
and without a MW 7r-pulse, Api, is plotted as a function of 
dark time r and fit to single-exponential decay. 



Fluorescence was collected by the same lens, passed 
through a dichroic mirror and was detected with a Si 
avalanche photodetector (Thorlabs APDllOA) in the 
range of 650-800 nm. A microwave (MW) field was ap- 
plied to the NV ensemble by means of a copper wire, 75 
/im in diameter, placed inside the cryostat in the vicin- 
ity of the optical focus. Light pulses were generated by 
passing cw laser light through two acousto-optical mod- 
ulators, which provided a combined extinction ratio of 
about 60 dB. MW pulses were produced by two series- 
connected switches, providing a combined isolation of 
120 dB. High isolation is critical, particularly for mea- 
surements of very long Ti , in order to avoid unintended 
changes in polarization from stray fields. Sequences of 
pulses were produced by a programmable pulse genera- 
tor PulseBlasterESR-PRO. 

Optical pumping and detection of spin polarization 
were performed using 532 nm excitation. Optical pump- 
ing occurs due to a spin-dependent intersystem cross- 
ing which transfers NV centers with spin projection, 
nis = ±1, to nis = ground-state sublevel [Fig. HJa)]. 
Due to the same process, NV centers in = 0, fluo- 
resce more brightly than NV centers in nis = ±1, al- 
lowing us to use the detected fluorescence intensity to 
determine the ground-state spin polarization. 

Figure [Ijb) shows the timing sequence used in the 
experiments to measure Ti. Each sequence begins with 
a 1 ms light pulse which polarizes NV centers into iris — 
0. We then apply a MW 7r-pulse, with typical duration 
on the order of 100 ns, which transfers NV centers into 
nis = 1 or TUs = — 1. Following a variable time delay, r, 
we apply another light pulse and detect NV fluorescence 
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FIG. 2. (a) 1/Ti as a function of temperature for S2, S3, 
and S8 samples. Symbols represent experimental data while 
solid lines are global fits to the model in Eq. ((T]). Error 
bars represent the standard errors from single-exponential 
fits. For consistency, we also plot previously reported exper- 
imental data [23,[2i|. (b) Ai [Eq. 0[ versus [NV[. (c) I/Ti 
(squares) and fluorescence intensity (circles) at 30 K versus 
position on the sample with a spatially- varying electron irra- 
diation dose. A magnetic field, B ~ 25 G was applied along 
[I00[ direction for all measurements presented in this figure. 



to determine the residual spin polarization. Afterwards 
we apply a control sequence with the 7r-pulse omitted, 
and subtract the two results. This technique results in 
a signal, Af;(t), which is proportional to the residual 
spin polarization, after time r, of only those NV centers 
which were excited by the MW 7r-pulse. Figure 1(c) 
shows a typical plot of Api versus r at T = 300 K, 
along with a flt to single-exponential function, A pi (r) = 
Ae^'^^'^^ . We found that all measurements reported in 
this manuscript can be satisfactorily fitted with a single 
decay constant. 

Insight into the nature of the longitudinal spin relax- 
ation mechanisms can be gathered from its temperature 
dependence. Figure [2] shows the relaxation rate l/Ti as 
a function of temperature for the three samples in Ta- 
ble H] A magnetic field, B w 25 G was applied parallel 
to the [100[ crystallographic orientation in order to sep- 
arate TUs = — >■ TUs — 1 and nig = — > nis = ~1 
transitions. At temperatures T > 200 K, 1/Ti is the 
same for all samples within a factor of « 2. For ex- 
ample, at 300 K, we measure Ti = 2.9(1), 5.5(1), and 
6.0(1) ms for S2, S8, and S3, respectively. This sug- 
gests Ti in this range is largely governed by intrinsic 
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processes such as interaction with phonons [26| . In con- 
trast, at low temperatures, a strong sample dependence 
of Ti is observed. The relaxation rate 1/Ti rapidly de- 
creases with decreasing temperature and flattens out at 
a certain value which is different for each sample. 

Considering interaction with lattice phonons, we note 
that the NV ground-state spin splitting, D « 3 GHz. 
Consequently for T ^ D/k w 0.1 K, where k is the 
Boltzmann constant, the occupation of phonon modes 
near 3 GHz is extremely low, and we can neglect relax- 
ation from a single-phonon process. However, higher- 
energy lattice phonons can produce Ti relaxation via 
a two-phonon Raman process. We expect this process 
to have the form (l/Ti) Raman oc [i^- Similarly, 
local phonons can produce Ti relaxation by resonant 
interaction with an excited vibrational level via a two- 
phonon Orbach-type process In this case, we ex- 
pect (l/Ti)orbach oc (e"feT — where A is the domi- 
nant local vibrational energy. 

Based on Fig. EJa), we also observe a leveling-off 
of (1/Ti) at low temperatures to a constant limiting 
rate. This rate is sample dependent and represents, as 
it will be discussed later, cross relaxation which arises 
from dipole-dipole interaction. We fit our data to the 
following global function: 
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where A2, A3, and A are fit parameters common to all 
samples, and ^i(S) is different for each sample. The fit 
is in excellent agreement with experimental data (see 
[2^ for a comparison with the fit parameters individual 
for each sample). The values for the fitting parameters 
are: Ai(S2) = 21(3) s"!, Ai(S8) = 0.06(1) s-\ Ai{S3) 
= 0.007(4) s-\ A2 = 2.1(6) X 103 s-^, A = 73(4) 
meV, and A3 = 2.2(5) x 10"" K-^g-i. The value of 
the parameter A closely matches the energy of the local 
vibrational modes of the NV center (28l . [29| . 

Figure El^b) shows ^1 as a function of NV concen- 
tration. The observed dependence is approximately lin- 
ear, and the fitted slope is 0.8(2) s~^ppm~^. To verify 
this observation we performed Ti measurements on a 
HPHT sample which has a spatial gradient in electron- 
irradiation dose and, correspondingly, in NV concentra- 
tion. The nitrogen concentration of this diamond was 
measured as in Ref. and found to be on the order 
of 100 ppm throughout the sample. We measured Ti at 
T = 30 K in four different spots along the dose gradi- 
ent. The relaxation rate 1/Ti and fluorescence intensity 
for each position on the diamond are presented in Fig. 
[2]Jc). Making an assumption that fluorescence inten- 
sity is proportional to NV concentration the results are 
consistent with linear dependence of Ai on [NV] . 

In order to gain further insight into the apparent cross 
relaxation processes, we measured magnetic fleld depen- 
dence of Ti for the sample with the highest impurity 
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FIG. 3. (a) l/Ti as a function of the magnetic field B applied 
along [111] direction (the corresponding transition frequen- 
cies are shown on the upper axis) for S2 sample at two differ- 
ent temperatures 296 K and 20 K. 1 /Ti was measured for the 
subensemble of NV centers aligned with B. Solid line rep- 
resents a fit (see text), (b) Transition frequencies between 
magnetic sublevels of NV centers and between sublevels of 
N as a function of magnetic field along ]111] direction. NV- 
NV cross relaxation occurs at G and 595 G, NV-N - at 514 
G. Solid lines correspond to the centers aligned with [111]; 
dashed lines correspond to the NV centers aligned at an ~ 
109° to this axis. 



concentration, S2. Magnetic fleld B was applied along 
the [111] axis using a neodymium ring magnet mounted 
on a 3-axis translation stage. For this magnetic fleld di- 
rection, one of four possible NV orientations is aligned 
parallel to B, while the other three orientations have 
symmetry axes which are at « 109° with respect to B. 
An additional 3-axis magnetic-coil system was used for 
flne tuning and scanning of the magnetic fleld. 

Figure|3](a) shows 1/Ti as a function of magnetic field 
for the subensemble of NV centers aligned along the 
magnetic field. The upper axis shows the corresponding 
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MW transition frequencies that were probed. The data 
were taken at two different temperatures: 20 K and 296 
K. Significant changes in 1/Ti are observed at certain 
magnetic field values: 0, 514, and 595 G. The relative 
changes are more pronounced for low-temperature data, 
where 1/Ti varies by more than an order of magnitude. 

Cross relaxation occurs when the transition frequency 
of two spin species, each having different spin polariza- 
tions, coincide. Transition frequencies as a function of B 
of the two sub-ensembles of NV centers (those co-linear 
and non-colinear with B) as well as that of the spin-1/2 
nitrogen impurities are shown in Fig. [3] (b). The 514 G 
peak apparently arises due to cross relaxation between 
NV centers co-aligned with B and the unpolarized ni- 
trogen spins. At and 595 G, the transition frequencies 
between the two subensembles of NV centers coincide. 
Cross relaxation can then occur if the subensembles are 
polarized differently by the initial optical pulse. For 
example, the peak at 595 G (1.2 GHz) occurs due to 
cross relaxation between jm^ = 0) — ^ |ms = — 1) tran- 
sition of NV centers aligned along the field direction 
and transitions between the eigenstates of other three 
degenerate orientations. At this field, each of the eigen- 
states is a mixture of the initial magnetic sublevels due 
to the large transverse field component. Cross relax- 
ation at G may be due to spin diffusion involving 
NV centers outside of the optical focus. Similar cross 
relaxation phenomena were previously studied for NV 
centers in diamond using other techniques, for exam- 
ple, monitoring changes of emission intensity [30- 32], 
measuring changes of the optical hole depth in the zero- 
phonon line [SS^ , and studying decay rate of stimulated 
spin-echo and spin- locking signals [30]. 

We fit our experimental results to three Lorentzians 
with central frequencies fixed at 0, 514, and 595 G [34"]. 
The amplitude at 20 K is a factor of two smaller than 
at 296 K for and 514 G peaks, and a factor of five 
smaller for the 595 G peak. We also observe a factor of 
two decrease in width for the 514 and 595 G peaks at 20 
K. We did not expect to see a temperature dependence 
of either amplitudes or widths (since dipolar coupling 
strengths are not expected to be temperature dependent 
in this range [l^); so the observed difference is a subject 
for future work. At both temperatures, the peak at 
514 G is the largest peak, which is consistent with the 
[N]:[NV] ratio in Table [B 

Let us look at the results from the point of view of 
magnetic sensing. Assuming that, ultimately, trans- 
verse relaxation times T2 on the order of Ti can be 
achieved [H-ill, (which may require the use of diamond 
samples depleted of ^'^C ^] and dynamic decoupling 
from the spin bath), the product [NVjTi is the figure of 
merit in terms of the sensor sensitivity Examining 
the data in Fig. Ufa), we find that at the lowest temper- 
atures, ]NV]Ti is roughly independent of concentration, 
owing to the aforementioned NV-NV cross relaxation. 
At higher temperatures, the relaxation rate is domi- 



nated by intrinsic phonon interactions and consequently 
independent of ]NV]. Thus the product ]NV]Ti increases 
linearly with concentration at a given temperature in 
this range. This motivates the use of high-density sam- 
ples in room-temperature applications. It is important 
to notice that, with appropriate choice of magnetic field, 
it is possible to "tune" the NV centers to a point with an 
optimum value of Ti. For example, applying bias field 
to a low-temperature high-density sample, it is possible 
to prolong Ti by more than an order of magnitude com- 
pared to near-zero field (Fig. [3] (a)). This property may 
also benefit other applications , including spin-ensemble 
memories in hybrid devices |36l439l |. 

In conclusion, we have investigated the temperature 
and magnetic-field dependence of longitudinal spin re- 
laxation of NV ensembles in diamond. We observed 
two distinct regimes (low and high temperature) , where 
different sets of relaxation mechanisms dominate. The 
dominant relaxation mechanism at high temperature is 
consistent with the exponential behavior of the two- 
phonon Orbach-type process and T^ behavior of the 
two-phonon Raman process. At low temperatures Ti is 
temperature independent and the relaxation is mainly 
determined by cross relaxation with neighboring spins. 
Through our magnetic-field studies, we identify this 
cross relaxation as being between differently polarized 
NV centers as well as between NV centers and nitrogen 
impurities. The longest Ti observed was 199(41) s at 
10 K for the CVD sample with low NV concentration. 
Our studies inform on the fundamental quantum limits 
of NV based computing and sensing applications. 
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TEMPERATURE-DEPENDENCE FITS FOR 
INDIVIDUAL SAMPLES 



In the main text, we presented a global fit of the 
temperature dependence of 1/Ti for the three samples 
investigated in this work. The global fit (Fig. 2(a) in the 
main text and Fig. ISlTa) here) describes well the 1/Ti 
variation over some five orders of magnitude. However, 
at the next ("factor-of-two") level of detail, the experi- 
mentally significant deviations for the high-temperature 
part of the dependence for the highest-NV-density S2 
sample, clearly seen in the insets of Fig. IS1[ become 
important. In order to understand these deviations, we 
performed fits of the temperature dependence given by 
Eq. (1) in the main text individually for each of the 
three samples studied in this work (Fig. ISlT b)). The 
values of the fit parameters are collected in Table ISll 
The fit results indicate that [within the model of Eq. (1)] 
there is an additional contribution to the relaxation ap- 
pears to be present in the highest-NV-density sample 
S2 at temperatures < 100 K. 

A detailed investigation into the nature of this ad- 
ditional relaxation awaits further work; however, as an 
initial step in that direction, we have measured 1/Ti 
at 296 K for different locations on the gradient sam- 
ple mentioned in the main text at magnetic field, B 
« 25 G applied parallel to the [100] crystallographic 
orientation (see Fig. 2(c)). The results are shown in 
Fig. [S2l The highest-NV-density part of the sample 
is similar in terms of its properties to the S2 sample. 
There appears to be a contribution to 1 /Ti that is cor- 
related with NV density, which one would not expect 
from Eq. (1) under the assumption that the coefficients 
in the second and third term are intrinsic to diamond. 
We note that the magnetic-field dependence of relax- 
ation shown in Fig. 3(a) in the main text also points 
towards a cross-relaxation contribution at room tem- 
perature exceeding what one would have expected from 
our simplified global picture. 



Ti COMPARISON WITH OTHER SYSTEMS 

In view of the importance of long-lived electron-spin 
systems in applications (which, in addition to mag- 
netometry, also include quantum information science. 



TABLE SI. Fit parameter values from fits of the Eq. (1) to 
the longitudinal spin relaxation data in Fig. Ell 



Sample 


Ai is-') 


A2 (s-^) 


A (meV) 


As (K'^s- 




Global fit 


S2 


21(3) 










S8 


0.06(1) 


2.1(6) X 10"^ 


73(4) 


2.2(5) X 10 


-11 


S3 


0.007(4) 










Individual fits 


S2 


20(3) 


4.7(2.4) X 10"'' 


77(6) 


2.1(8) X 10 


-11 


S8 


0.06(1) 


1.7(6) X 10"^ 


69(4) 


2.5(5) X 10 


-11 


S3 


0.007(3) 


2.1(4) X 10"^ 


76(4) 


2.2(3) X 10 


-11 



where the electron-spin Ti is a limiting factor P, 0[ , in 
Table IS2| we list representative values of Ti for both 
cryogenic and room temperatures. The Table clearly 
illustrates the advantages of the nitrogen vacancy cen- 
ters, especially, at room temperature. 



TABLE S2. Comparison of approximate longitudinal spin 
relaxation time values for different defects in solids which 
are of high interest for spin-based information processing 
devices 



Defect type 



T (K) B (G) Ti (s) Reference 



NV in diamond 300 

NV in diamond 10 

Defects in 4H-SiC 300 

Ga(In)As quantum dot 1 

Si/Si02 quantum dot 0.05 

Si/SiGe quantum dot 0.05 

P donors in Si 10 

P donors in Si 1.8 



25 
25 


4 X lO"* 
2 X lO"* 

1.5 X 10'' 

3.5 X 

3.5 X 



6 X 10" 



lO'' 
10^ 



2 X 10^ 

2 X 10"'' 

2 
2 
1 

3 

|3 



X 

2 X 10" 
4 X 10" 
2 X 10" 

1 X 10" 

2 X 10^ 



This work 
This work 
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m 
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FIG. SI. 1/Ti as a function of temperature for S2, S3, and 
S8 samples. Symbols represent experimental data while solid 
lines are fits to the model in Eq. (1) where (a) all param- 
eters except Ai(S) are common for all samples and (b) all 
parameters are variable. 
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FIG. S2. (a) l/Ti and (b) fluorescence intensity at 296 K 
versus position on the sample with a spatially- varying elec- 
tron irradiation dose. 



